A Lattice Test of 1/iVc Baryon Mass Relations 



o 
o 

(N 

(D 

Q 
o 



in 
I 

Oh' 



> 
On 
(N 

vn 
o 

l> 
o 

ON 

o 



- '—I 



Elizabeth E. Jenkins.^ Anccsh V. Manohar,^ J.W. Negele,^ and Andre Walker-Loud'^ 

^Department of Physics, University of California at San Diego, La Jolla, CA 92093 
^Center for Theoretical Physics, Massachusetts Institute of Technology, Cambridge, MA 02139 
^Department of Physics, College of William and Mary, Williamsburg, VA 23187-8795 
(Dated: December 30, 2009 21:14) 

1 /Nc baryon mass relations are compared with lattice simulations of baryon masses using diflerent 
values of the light-quark masses, and hence different values of 517(3) flavor-symmetry breaking. 
The lattice data clearly display both the 1/Nc and SU{3) flavor-symmetry breaking hierarchies. 
The validity of 1/Nc baryon mass relations derived without assuming approximate SU{3) flavor- 
symmetry also can be tested by lattice data at very large values of the strange quark mass. The 1/Nc 
expansion constrains the form of discretization effects; these are suppressed by powers of l/N^ by 
taking suitable combinations of masses. This 1/Nc scaling is explicitly demonstrated in the present 
work. 



I. INTRODUCTION 

The 1/Nc expansion of QCD [1] is a valuable tool for 
studying the nonperturbative dynamics of the strong in- 
teractions [2, 3]. In the limit oo, the baryon 
sector of QCD has an exact contracted SU{2Np) spin- 
flavor-symmetry [4]. For finite N^, the contracted spin- 
flavor-symmetry is broken by effects suppressed by pow- 
ers of 1/A''c [4, 5]. The spin-flavor structure of the 1/Nc 
breaking terms is predicted at each order in the 1 /Nc ex- 
pansion [4-6] . The spin-flavor structure of many baryon 
properties have been derived in a systematic expansion 
in 1/Nc [7-12], and the results are in excellent agreement 
with experiment (for reviews, see [13, 14]). 

One important application of the baryon 1/Nc expan- 
sion is to baryon masses [6, 15, 16]. By choosing ap- 
propriate linear combinations of the baryon masses, one 
can study coefRcients of the baryon mass 1/Nc opera- 
tor expansion with definite spin and flavor transforma- 
tion properties. In the case of perturbativc SU (3) flavor- 
symmetry breaking, the 1/Nc analysis gives a hierarchy 
of baryon mass relations in powers of 1/Nc and the di- 
mensionless SU{3) breaking parameter e oc rUs/A^ [15]. 
The analysis in Ref. [15] showed that the experimentally 
measured masses of the ground state octet and decuplet 
baryons exhibit the predicted hierarchy of the combined 
expansion in 1/Nc and SU{3) flavor-symmetry breaking. 
The 1/Nc expansion also has been used to obtain very 
accurate predictions for the charm and bottom baryon 
masses [16] (to better than 10 MeV accuracy) which 
have been confirmed by recent experiments. In addition, 
Ref. [6] derived baryon mass relations which only depend 
on the 1 /Nc expansion and which arc valid even if fiavor 
SU{5) is not a good approximate symmetry, i.e. for large 
values of the strange quark mass. 

The predictions of the 1/Nc expansion for baryon 
masses are in excellent agreement with the experimen- 
tal values. However, in comparing with the experimental 
data, one is restricted to only one value of Nc and to the 
physical quark mass values. Testing the predictions as 
a function of light-quark masses and Nc is now possible 



with very accurate simulations of baryons using lattice 
QCD. 

Tremendous progress in lattice QCD has been achieved 
recently in the simulation of baryon masses using differ- 
ent values of the light-quark masses. Extrapolation of 
baryon masses on the lattice to the physical point has re- 
produced the experimental values at the 1 — 3% level [17]. 
The lattice data, however, contain important additional 
information about the dependence of the baryon masses 
on the quark masses, which can be utilized. Simula- 
tions of baryon masses have been performed as a function 
of SU (3) flavor-symmetry breaking ranging from small 
perturbative flavor-symmetry breaking to large nonper- 
turbative flavor-symmetry breaking. There arc also lat- 
tice simulations at different values of Nc (though not for 
baryons) which are able to test Nc scaling rules [19, 20]. 

In this paper, we show that existing lattice simulations 
(at Nc = 3) of the ground state baryon masses already 
are sufficiently accurate to exhibit and test interesting 
features of the 1/Nc and SU (3) flavor-symmetry breaking 
expansions. Still more accurate simulations are needed 
to test the most suppressed mass combinations of the 
1/Nc expansion, but continued improvements in lattice 
simulations of baryon masses are expected in the short 
and long term, so it should eventually be possible to test 
these relations as well. We discuss in this paper how 
present and future lattice data can be utilized to study 
the spin and flavor structure of baryon masses. Although 
we do not focus on this point here, it should eventually be 
possible to test the Nc scaling rules in the baryon sector 
by lattice simulations which vary the number of colors 
Nc- 

An important observation is that the 1/Nc counting 
rules hold at finite lattice spacing, and so are respected 
by the lattice results including the finite lattice spacing 
corrections dependent on the lattice spacing a. Thus the 
discretization corrections are constrained by the 1 /Nc ex- 
pansion. 

Lattice computations of hadron masses are done at 
varying values of the light-quark masses m„, rud and m^, 
usually in the isospin limit to„ = = rriud- The lat- 
tice results are then extrapolated to the physical values 
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FIG. 1: One-loop correction to the baryon mass due to ■r,K 
radiative corrections. 



of TTiud- It has not been possible to compute hadron 
masses for physical values of m„(; yet due to the large 
computational time needed, since Tcomp oc l/m^^.^ How- 
ever, recently, with both algorithmic developments [21] 
and large parallel computing machines, there are two 
groups simulating at or near the physical light-quark 
mass point [17, 22]. 

The light-quark mass dependence of the hadron masses 
is determined by chiral perturbation theory. There are 
nonanalytic in contributions from loop corrections 
which are calculable, as well as analytic terms which de- 
pend on low-energy constants (LECs) of the chiral La- 
grangian. In the baryon sector, the leading nonanalytic 
terms are m'J'^ and m^lumq. The m^^^ term, which 
arises from Fig. 1 and is proportional to / {Wn f'^) 

times axial couplings, is found to be rather large, naively 
of order a few hundred MeV. This term is absent for the 
pseudo-Goldstone boson masses, but is present for other 
mesons such as the vector mesons [23]. 

A surprising feature of recent lattice results is that the 
baryon masses as a function of MTr,K,ri do not show a 
large dependence. Fits to the baryon masses also 
give much smaller values for the baryon axial couplings, 

Q 3/2 

and require that the AI ~ niq term is almost com- 
pletely cancelled by Inm^ and terms. This cancel- 
lation must be accidental at the currently accessible lat- 
tice quark masses, since the terms have different rriq de- 
pendence. An alternative conclusion is that lattice quark 
masses are too large for 5*^7(3) chiral perturbation the- 
ory to be valid, and that perturbative chiral behavior sets 
in only for much smaller quark masses than the strange 
quark mass. This conclusion, however, fails to explain 
why SU{3) flavor-symmetry is so evident in baryon phe- 
nomenology. 

The 1 /Nc expansion constrains the structure of baryon 
chiral perturbation theory [4, 6, 9, 18]. Chiral corrections 
have to respect the spin-flavor structure of the ex- 
pansion, and this leads to some important restrictions 
on the form of the chiral loop corrections. For exam- 



One inverse power of the quark mass is from the conjugate- 
gradient inversion of the fermion Dirac operator, which scales 
with the condition number. The additional inverse powers of 
the light-quark mass arise from estimations of (i) the increased 
auto-correlation time of the HMC evolution as well as (m) the 
molecular dynamics step-size used to adjust the acceptance rate. 
There are additional costs not counted here. For example, to 
keep the volume corrections exponentially small, one must work 
with rriTrL > 4, where L is the spatial size of the lattice. 



pie, the baryon mass, which is of order Nc [2, 3], gets 

3/2 

nonanalytic corrections proportional to TOs , which are 
large. This large nonanalytic contribution might lead one 
to expect that there should be large deviations from the 
Gell-Mann-Okubo mass relations, which were derived as- 
suming that the mass operator was linear in m^. One 
can show, however, that the NciTis^^ terms are a spin- 
flavor singlet, and give the same contribution to octet 
and decuplet baryons, whereas the m^^"^ terms are a 

spin-singlet flavor-octet, and the mV^ /Nc terms are spin- 
singlet flavor-27. Only the latter terms contribute to the 
Gell-Mann-Okubo mass combinations, so that deviations 
from these relations are a factor of l/N"^ smaller than 
naive expectation and consistent with experiment. The 
small size of the Gell-Mann-Okubo relation was recently 
confirmed for a range of light-quark masses [24]. For 
QCD, the cancellation to order 1/Nc arises as a numeri- 
cal cancellation between octet and decuplet intermediate 
states [25] ; to see the parametric form of the cancellation 
in \/Nc requires computing the chiral corrections using 
the SU (3) flavor representations of baryons containing 
Nc quarks. 

The mass relations of Ref. [15] project the baryon 
masses onto different spin-flavor channels. By studying 
these mass relations as a function of m^, one can in- 
vestigate whether unexpected chiral behavior arises in a 
particular channel. 

The organization of this paper is as follows: In Sec. II, 
the baryon mass relations of the 1 /Nc expansion are sum- 
marized briefly, and in Sec. Ill, the lattice simulation data 
are described. Sec. IV presents the results of a computa- 
tion of 1/Nc mass combinations on the lattice for vary- 
ing values of S'[/(3) flavor-symmetry breaking. Sec. V 
discusses the lattice analysis using heavy baryon chiral 
perturbation theory [25]. Our conclusions are presented 
in Sec. VI. 



II. 1/N^ BARYON MASS RELATIONS 

The \/Nc expansion of the baryon mass operator for 
perturbative S'[/(3) flavor-symmetry breaking is^ 

AI - A/1-0 + A/«-0 + Af27.o + M^^^\ 
AI^^ = c^Nct + c\^j^^ 

M''''^c%fl-^{TMT\T']], (1) 



^ We use the notation and conventions of Ref. [15] . 
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where the superscript denotes the flavor SU{S) represen- 
tation and the spin SU (2) representation of each term. 
The 1/Nc expansion of the baryon mass operator con- 
tains eight independent operators, corresponding to the 
eight isospin multiplets of the ground state baryons. The 
arbitrary coefficients C(„) multiplying the l/Nc operators 
are functions of 1 /Nc and SU (3) flavor-symmetry break- 
ing. Each coefhcicnt is order 1 at leading order in the 
l/Nc expansion. The non-trivial content of the l/Nc ex- 
pansion is the l/Nc suppression factors for the differ- 
ent terms in Eq. (1). Following Ref. [15], we define the 
parameter e oc ttIs/A^, which is a dimensionless mea- 
sure of flavor-S'?7(3) symmetry breaking. SU{3) flavor- 
symmetry breaking transforms as a flavor-octet, so the 1, 
8, 27 and 64 operator expansions are zeroth, first, sec- 
ond, and third order in SU (3) flavor-symmetry breaking, 
respectively. Thus, the coefficients c^'", c^^'" and c^"''" 
are naively order e, and e^, respectively, in the SU{3) 
flavor-symmetry breaking expansion. 

The 1/A^c expansion for the baryon mass operator for 
SU{2) X [/(I) flavor-symmetry [6] is given by 



M = Ncl 
1 



K - 
1 



1 

'n'c 



-I' 



Nc 



■J' 



Nc 



N? 



(2) 



where each operator is accompanied by an unknown co- 
efficient which has been suppressed for simplicity. Each 
coefficient is order unity at leading order in the 1/A^c 
expansion. The operator K = Ns/2, where Ng is the 
s-quark number operator. Again, there are eight inde- 
pendent operators in the l/Nc expansion corresponding 
to the masses of the eight isomultiplcts of the ground 
state baryons. 

The l/Nc expansion continues to hold at flnitc lat- 
tice spacing, so that expansions of the form Eqs. (1) and 
(2) are valid with coefficients that depend on the lattice 
spacing a.^ Thus, the order Nc discretization correction 
is universal, and has the same value for all the octet and 
decuplct baryons. This term can be eliminated by study- 
ing baryon mass ratios or differences. The largest term 
that produces different discretization effects for the octet 
baryons is the term in M^'", or equivalently the K 
term in Eq. (2), which results in an 0(1) mass splitting 
that is suppressed by one factor of l/Nc relative to the 
leading mass term Net. 

The twelve 1 /Nc mass relations that we study are tab- 
ulated in Table I. The eight mass relations Mi-Mg corre- 
spond to the first eight mass combinations in Table II of 
Ref. [15] . These mass relations are the mass combinations 
which occur at definite orders in l/Nc and perturbative 
SU{3) flavor-symmetry breaking; each relation picks out 



^ Rotational symmetry is broken down to a discrete cubic symme- 
try, so that irreducible representations of the spin 5(7(2) group 
are replaced by irreducible representations of the cubic group. 



a particular coefficient in Eq. (1). The mass combina- 
tions in Table I correspond to the coefficients listed in 
the table. There are two mass combinations, Mi and 
M2, which are SU{3) singlets and which occur at orders 
Nc and 1/A^c, respectively. There are three flavor-octet 
mass combinations, M3, M4 and M5, which arc propor- 
tional to one factor of SU (3) symmetry breaking e and 
which occur at orders 1, l/Nc and l/N^, respectively. 
There are two flavor- 2 7 mass combinations, Mg and M7, 
which occur at second order in SU{3) breaking and at 
orders l/Nc and 1/N^ in the l/Nc expansion. Finally, 
there is one mass combination Mg which is suppressed 
by three powers of SU{3) breaking and by l/Nc 
the 1/A^c expansion. The additional four mass combina- 
tions Ma^M d correspond to the flrst four mass relations 
of Ref. [6]. These baryon mass combinations arc each 
order l/Nc the l/Nc expansion. Only three of the 
four combinations Ma-Mo arc linearly independent, cor- 
responding to the fact that there are three independent 
operators at order 1/N^ in Eq. (2). These mass relations 
were derived assuming only isospin flavor-symmetry, and 
so do not assume approximate SU (3) flavor-symmetry 
and are valid even for very large SU(3) flavor-symmetry 
breaking. Note that the l/Nc mass relations M^-M^i 
are related to the l/Nc relations M5, M7 and Mg by 



ivu 
Mb 
Mc 



ID = 



1 1 2 

— M5 H M7 Mg, 

10 70 7 

14 2 

M.; H M7 Mg, 

15 105 21 ^' 
-^(M7 + Mg), 

--Mg. 

2 ^ 



(3) 



Each mass relation deflnes a mass combination of the 
octet and decuplet masses Mi of the form 



M = ^ c,Mi 



(4) 



In this work, isospin breaking is neglected, so each Mi 
denotes the average mass of a baryon isomultiplet. The 
normalization of each mass relation is arbitrary, and de- 
pends on a particular choice of normalization for coeffi- 
cients in the Hamiltonian, i.e. c, Xci defines another 
mass relation with the same spin and flavor quantum 
numbers. To remove the rescaling ambiguity, Ref. [15] 
used the accuracy defined by 



A 



(5) 



to quantify the fractional error of a given mass relation. 
Here we use a related quantity, the scale invariant mass 
combination 



R 



(6) 
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TABLE I: Mass combinations Mi-Mg from Ref. [15] and Ma-Md from Ref. [6]. The coefficients and orders in 1/Nc and 
perturbative 5f/(3) flavor symmetry breaking e are given for mass combinations Mi-Mg. Combinations Ma-M_d are obtained 
at order l/N^ assuming only isospin flavor-symmetry. 



Label 


Operator 


CoefBcierit 


Mass Combination 


1/Nc 




Ml 

Ma 
Ms 

M4 

Ms 
Me 
Mt 
Ms 


1 

{J^^8} 

{T ,T } 

{r*,{r®,T**}} 


160 iV, cJq° 

1 oc 1 2 „27,0 
_63 1 £2 27,0 
a./3 1 ^3 64,0 


25(2Af + A + 3E + 2H) - 4(4A + 3E* + 2H* + Q) 
5{2N + A + 3E + 2H) - 4(4A + 3E* + 2H* + Q.) 
5{6N + A - 3E - 4H) - 2(2A - E* - Q) 

iV + A - 3E + H 
(-2iV + 3A - 9E + 8H) + 2(2A - S* - Q) 
35(2A — 3A — E + 2^j — 4(4A — 5E — 2^ + 3s 2) 
7(27V - 3A - E + 2H) - 2(4A - 5E* - 2H* + SQ) 
A - 3E* + 3H* - Q 


1 

1/iVc 

l/Nc 
1/m 


1 
1 

e 
e 
e 

2 

e 


Ua 
Ms 
Mc 
Md 






(E*-E)-(H*-H) 
i(E + 2E*)-A-|(A-7V) 
-i (2Af - 3A - E + 2H) + i (A - E* - H* + !^) 
-i (A-3E* +3H* - fi) 







Dividing by |ci| instead of by i |ci|A/i avoids mix- 
ing difFercnt flavor representations via the denominator 
factor. The rescaled relations Ri-Rg and Ra-Rd have 
dimensions of mass. 

In our numerical analysis, we shall use the dimension- 
less variable 

e = (7) 



respectively, on the m007, mOlO, m020, m030, m040, and 
m050 ensembles, where the labels denote the light-quark 
masses in lattice units. ^ In the dynamical ensembles and 
the computation of the valence propagators, the strange 
quark was held fixed near its physical value. (In fact the 
strange quark was ~ 25% too large [36].) For further 
details of the calculation, we refer the reader to Ref. [27]. 



as a measure of SU{i) breaking, where A^ ^ AttJ ~ 
1 GeV [26] is the scale of chiral symmetry breaking. 



III. LATTICE SIMULATION 

In this work, we use the results of the recent LHPC 
spectrum calculation [27] to explore the mass combina- 
tions of the l/Nc expansion. The LHP Collaboration 
utilized a mixed-action lattice calculation with domain- 
wall [28-30] valence propagators computed with the Asq- 
tad improved [31, 32] dynamical MILC gauge ensem- 
bles [33, 34]. The calculation was performed at one 
lattice spacing with a ^ 0.125 fm, and a fixed spatial 
volume L ^ 2.5 fm. The pion and kaon masses used 
in Ref. [27] are {M^,Mk} = {293,586}, {356,604}, 
{496,647}, {597,686}, {689,729} and {758,758} MeV, 



Using the bootstrap resampled lattice data, we deter- 
mine the 12 mass combinations of Table I on each en- 
semble. The results arc collected in Table II. These 
results were determined with the absolute scale of 
= 1588 MeV on all coarse ensembles, where the 
scale used in Ref. [27] was determined from heavy- 
quark spectroscopy. We have additionally determined 
the mass combinations using the smoothed values of 
ri/a, where ri is determined on each different ensemble 
from the heavy-quark potential with r^i^(ri) = 1 [37]. 
The values of determined in this way range from 
{1597, 1590, 1614, 1621, 1628, 1634} MeV from the light- 
est to heaviest quark mass. These two scale-setting meth- 
ods are in good agreement, as shown in the next section. 



* The strange quark and many of the Ught-quark propagators were ^ The mass Mr) is defined at this order by the Gell-Mann— Okubo 
computed by the NPLQCD Collaboration [35]. formula = |A/2, - ^M^. 
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TABLE II; Values of mass combinations Mi in GeV using the scale setting a ^ — 1.588 GeV. 





0.2579 


0.2380 


0.1718 


0.1141 


0.0574 





Ml [GeV] 


192(1) 


197(1) 


211(1) 


221(2) 


237(2) 


242(1) 


Ma [GeV] 


-12.2(5) 


-13.3(6) 


-10.9(4) 


-10.6(6) 


-7.8(6) 


-8.2(4) 


M3[GeV] 


-8.07(20) 


-7.15(15) 


-4.40(08) 


-2.81(09) 


-1.28(04) 





M4[GeV] 


-0.214(10) 


-0.181(7) 


-0.099(5) 


-0.056(6) 


-0.022(3) 





Ms [GeV] 


0.29(7) 


0.35(8) 


0.13(3) 


0.14(5) 


0.04(2) 





Ms [GeV] 


-1.05(36) 


-0.25(26) 


-0.003(71) 


0.15(14) 


-0.02(1) 





M7[GeV] 


-0.30(12) 


-0.05(11) 


0.005(24) 


0.08(7) 


-0.003(4) 





M8[GeV] 


0.02(1) 


0.012(09) 


-0.001(1) 


0.015(12) 


0.00004(6) 





Ma [GeV] 


0.0004(58) 


0.018(7) 


0.001(3) 


0.004(4) 


-0.0003(19) 





Ms [GeV] 


-0.022(10) 


-0.018(10) 


-0.0004(30) 


-0.003(3) 


0.0000(13) 





Mc[GeV] 


0.010(4) 


0.002(4) 


-0.0002(8) 


-0.003(3) 


0.0001(1) 





MD[GeV] 


-0.012(5) 


-0.0061(45) 


0.0004(6) 


-0.0076(58) 


-0.00002(3) 







0.2 0.4 0.6 

M„ (GeV) 



FIG. 2: Octet and decuplet baryon lattice masses as a func- 
tion of Mtt. The blue squares are the PDG values and the 
triangles are LHP Collaboration lattice data simulated with 
fixed rus and varying m„d- The two right-most points sim- 
ulate the flavor-singlet baryon masses in the exact 5(7(3) 
flavor-symmetry limit M-n = Mk- Note that the value 
of Mtt at the SU{3) flavor-symmetry point is quite large, 
= 758 MeV. The (PDG) baryons from bottom to top 
are A^, A, E, A, H, E*, H* and Q. 



IV. COMPARISON WITH LATTICE DATA 

A plot of the lattice baryon masses as a function of 
Mtt is given in Fig. 2. The triangles are the lattice data, 
and the blue squares are the Particle Data Group (PDG) 
values.^ The usual way of studying the data is to fit to 
the individual baryon masses. The analysis shows 



Note that Mk varies slightly for the different M-^ values and is 
not exactly equal to the physical M^- 



that it is useful to study the mass combinations in Ta- 
ble I instead, as these combinations have definite orders 
in a-nd SU{3) symmetry breaking. 

In this section, we present a series of plots of the 12 
scale invariant baryon mass combinations Ri-Rg and 
Ra-Rd as a function of SU{i) flavor-symmetry break- 
ing. In each plot, the red triangles and green circles are 
the lattice data, and the errors on the mass combinations 
have been computed using the bootstrap data sets of the 
octet and decuplet masses, taking advantage of the full 
correlations in the lattice data. The red triangles and 
green circles are the lattice results using the two differ- 
ent scale-setting methods. The red points are determined 
with the absolute scale setting = 1588 MeV, and the 
green points are determined with the smoothed ri /a val- 
ues [37]. The blue square is obtained using the physical 
baryon and meson masses from the PDG [38], with the 
error bar computed using experimental uncertainties on 
the isospin averaged masses. For most of the plots, the 
error on the blue point is smaller than size of the point, 
and is not visible. 

The horizontal axis is — in units of (GeV)^, 
which is a measure of S't/(3) flavor-symmetry breaking, 
and the vertical axis is the mass combination in MeV. 
For mass combinations proportional to powers of SU (3) 
breaking, we plot both the mass combination and the 
mass combination divided by the appropriate power of e, 
using the definition Eq. (7) for e, so that the units remain 
MeV. 

The average 0{Nc) mass of the octet and decuplet 
baryons is ^ 1000 MeV. Naive power counting with 
A^c — 3 implies order Nc masses are ~ 1000 MeV, order 1 
masses are ^ 300 MeV, order 1/iVc masses are ^ 100 MeV, 
and order l/N"^ masses are ~30 MeV. This 1/A^c hierar- 
chy is evident in the PDG and lattice data for the baryon 
mass combinations Ri-Rs and Ra-Rd below. 

Figure 3 plots the first relation which is the scale 
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1500 



1000 
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0.1 0.2 0.3 
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FIG. 3: Ri as a function of {M% ~ Ml). The red triangles 
and green circles are lattice results with two different scale 
setting methods (see text), and the blue square is using the 
PDG values for the hadron masses. R\ is the average 0{Nc) 
mass of the ground state baryons. 



FIG. 4: i?2 as a function of {Mj^ ~ Ml). The red triangles 
and green circles are lattice results with two different scale 
setting methods (see text), and the blue square is using the 
PDG values for the hadron masses. R2 is minus the 0{\/Nc) 
hyperfine mass splitting between the spin-3/2 decuplet and 
the spin-1/2 octet baryons. 



invariant mass corresponding to the iVcl operator in the 
1/Nc expansion. This mass relation is order 1000 McV 
over the range of SU{2i) flavor-symmetry breaking simu- 
lated by the LHP Collaboration and is seen to be fairly 
independent of 5L/(3) breaking as expected. The devia- 
tion of the physical PDG point from the lattice data is 
presumably due to the larger than physical strange quark 
mass of the simulation. In principle, this relation also is 
susceptible to finite discretization errors, beginning at 
0{a?), although these are expected to be small for this 
data set [27]. 

Figure 4 plots the second relation i?2, which corre- 
sponds to the J^/iVc operator in the 1/Nc expansion. 
The mass combination i?2 is clearly suppressed relative 
to the average ground state baryon mass This hy- 
perfine splitting is predicted to be order 1/Nc, and so of 
order 100 MeV in magnitude. The lattice data allow us to 
study this mass relation as a function of SU (3) breaking. 
Notice that the relation works independently of SU{3) 
breaking as predicted theoretically, and that R2 does not 
vanish at the SU{3) symmetry point A/|- = M^. The 
numerical results clearly show that the suppression of R2 
is a consequence of 1/iVc, and not due to a hidden flavor 
breaking suppression factor. 

Figure 5 plots the third relation i?3, which is of order e. 
The relation vanishes as M]^ — Ml — > 0, as can be seen 
in the upper panel. The lower panel divides the mass 
relation by e, and shows that R^/e is an 0(1) mass in 
the 1/Nc expansion, or ~ 300 MeV in magnitude. Notice 
that R^/e does not vanish in the 5'C/(3) symmetry limit. 
The SU (3) symmetry point of the lattice data is omitted 
in the R^/e plot, as it involves dividing by zero. 

Figure 6 plots i?4 and i?4/e, which are of order e/Nc 
and 1/Nc, respectively. The mass combinations i?3 and 
i?4 are both order e, but i?4 is suppressed relative to 
i?3 by an additional factor of l/Nc- This suppres- 



sion is clearly visible in the numerical values which are 
~ 100 MeV in magnitude. i?4 is seen to vanish in the 
SU{i) symmetry limit, whereas R^/e does not vanish at 
the SU{i) symmetric point. Both Rz/e and i?4/e display 
some dependence on S'C/(3) symmetry breaking, which 
implies that there is e dependence in the coefficients c^]^^ 

and c^2), respectively. 

Figure 7 plots R^ and R^/e, which are of order e/N^ 
and 1/Nl, respectively. The mass relation R^/e is con- 
sistent with being a 1 /iV^ mass of magnitude 30 MeV. 
The lattice data points now have relatively large errors, 
however. 

i?6, Ri and i?8 are plotted in Figs. 8, 9 and 10, respec- 
tively, for completeness. In each case, the physical point 
(blue square) has the expected size, and the lattice results 
arc compatible with 1 /Nc and e power counting expecta- 
tions, but the error bars are now rather large. It would 
be very interesting to have more precise lattice data for 
the Rq and i?7 mass relations because the leading SU (3) 
fiavor-symmetry breaking contribution to the flavor-27 
mass splittings in chiral perturbation theory is 0(e^/^) 
rather than the naive 0{e^) of second order perturbative 
fiavor-symmetry breaking. With very precise data, as in 
Ref. [39], it would be possible to test this prediction of 
chiral perturbation theory. 

The last four plots. Figs. 11-14, are the relations 
Ra-Rd, which are of order l/N^ and assume only 
SU{2) X C/(l) flavor-symmetry. From the plots, it is evi- 
dent that Ra-Rd are satisfied irrespective of the value of 
SU{3) breaking. Ra-Rd are aU predicted to be 0{1/Nl) 
masses, or approximately 30 MeV. The figures show that 
Ra-Rd are significantly smaller than 30 MeV. This sup- 
pression can be understood from the perturbative SU (3) 
fiavor-symmetry breaking analysis, since the lattice data 
varies e only over the range [0, 0.2579] and does not con- 
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FIG. 5: Ri and Rs/e as a function of {Mj^ - M^). Rs is 
0{e) at leading order in the 1/Nc and SU{3) flavor-symmetry 
breaking expansions, and Ra/e is 0(1). 



FIG. 7: Rb and Rs/e as a function of (Aljc - Ml). Rs is 
0{e/Nl) at leading order, and R^/e is ©(l/iV^^). 
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M,' - M; (GeV^) 



FIG. 6: 7?4 and Ri/e as a function of {M]( - Ml). Ra is 
0{t/Nc) at leading order, and Ri/e is 0{1/Nc). 



FIG. 8: Ra and iJs/e^ as a function of (A/|- - Af^). Tie is 
©(eViVc), and iie/e' is 0{l/Nc). 
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FIG. 9; i?7 and Et/e^ as a function of {M'k - Ml). Rr is 
0{e^/N^), and Rr/e'^ is 0{1/Nl). 




0.3 



M,' - M; (GeV^) 



FIG. 10: Rs and iJs/e as a function of {M% - Ml). Re. is 
O(eVA^c), and R^/e' is 0(l/iV,^). 




M; - M; (GeV^) 



FIG. 11: Ra as a function of {M'i - M^). 7?a is 0{l/Nl 




M; - M; (GeV^) 



FIG. 12: iifl as a function of {M^i - Ml). Rb is 0{1/Nl). 



tain large nonperturbative values of e. From Eq. (3), one 
sees that for small values of e, Ma and Mb are 0{e/N^) 
since they are linear combinations containing M5; Mc is 
0{e^ /N^) since it is a linear combination containing M7; 
and Md is 0{e^/N^) since it is proportional to Mg. It 
would be interesting to see how well relations Ra-Rd 
work at even larger values of SU (3) breaking extending 
to very large values for which the SU (3) flavor-symmetry 
is completely broken. In the limit — > 00, relation Ra 
goes over to the mass relation for heavy-quark baryons, 
Sg — Eq = Sq — Sq, so it would be interesting to look 
at the crossover from light-quark mass behavior propor- 
tional to rris to heavy-quark behavior in which the hyper- 
fine splitting is proportional to l/mg. Similarly, in the 
limit rus — > 00, relation Rb goes over to the heavy-quark 
baryon mass relation [i (Sg 4- 2Sg) - Ag] -| (A - iV), 
which is independent of the heavy-quark mass. The 
crossover from ~ rrig behavior to mass-independence in 
the rris —>■ 00 limit also would be interesting to observe. 
In addition, it would be interesting to determine the value 
of e at which the perturbative SU{3) flavor-symmetry 
breaking analysis fails. Finally, an extended range for 
e oc M|- — Ml would allow one to see whether relations 
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M; - M; (GeV^) 



FIG. 13: Rc as a function of {Ml -Ml). Rc is 0{1/N'i 




M; - M; (GeV^) 



FIG. 14: Rd as a function of {M^ - Ml). Rd is 0{l/N^ 



Ra-Rd remain independent of e even for large values of 
e. 



V. HEAVY BARYON xPT ANALYSIS 

As mentioned in the introduction, a recent surprise 
from current lattice light-quark spectrum calculations is 
that the baryon masses do not show a large ^ ^ be- 
havior. In fact, the baryon masses display an unexpected 
(nearly) linear in scaling for a large range of pion 
masses [27]. In the case of the nucleon, this finding was 
verified in all current 2-1-1 flavor lattice calculations of 
the nucleon mass [40]. These findings pose an interest- 
ing theoretical puzzle and indicate the presence of large 
cancellations to the baryon spectrum from the different 
orders in the chiral expansion. Further, recent SU (3) chi- 
ral extrapolations of the baryon masses [40, 41], carried 
out using the next-to-leading order (NLO) extrapolation 
formulae, have found that the axial couplings, D, F, C 
and H, when left as free parameters, are returned from 
the analysis with values significantly different from those 
expected based upon phenomenology [42] or the recent 



lattice determination [43]. It is not clear at this point 
whether these findings are due to a breakdown of SU (3) 
heavy baryon xPT in this mass range, or, for example, 
the necessity of including the next-to-next-to-leading or- 
der (NNLO) terms in the analysis. 

In this section, we perform an SU{S) heavy baryon 
xPT extrapolation analysis of the various mass relations 
Mj presented in Table II. This allows us both to extrap- 
olate our results to the physical {Af^,Af/f} limit and 
to begin exploring the combined l/Nc and SU{3) ex- 
pansion utilizing precision lattice data. An important 
question to explore is whether this combined expansion 
displays significantly improved convergence properties to 
the SU{3) extrapolations performed previously. Addi- 
tionally, because the different mass relations are sensitive 
to the different coefficients of the baryon mass expansion 
in Eq. (1), we can directly extrapolate different orders in 
the SU (3) chiral expansion. For example, as mentioned 
above, the leading contribution to Mg and M7 scales as 
0(6"^/^) in heavy baryon xPT. Finding a definitive sig- 
nal of this scaling would be a significant confirmation of 
the nonanalytic light-quark mass behavior predicted by 
chiral perturbation theory. 

The light and strange quark mass dependence is de- 
scribed by SU{5) heavy baryon chiral perturbation the- 
ory (HBxPT) [25] . The baryon chiral Lagrangian is given 
by 

jC = Tr {Biv ■ VB) - [iv • P - Aq] 

+ 6z5Tr {B{M+,B}) + bpTr {B[M+,B]) 

+ boTr{BB)TT{M+) 

+ 7A/ f^M+T^ - a T'^T^Tr {M+) 

-f 2D Tr {BS^'iA^,, B}) + 2F Tr {BS^iA^, B]) 

+ 2n fS^A-'Tf, + C {ft'A^B + BA^.T") , (8) 

where B and are the octet and decuplet fields, re- 
spectively; Vf^ is the four- velocity of the baryon and 5^ is 
the spin- vector. The decuplet-octet mass splitting in the 
chiral limit is Aq, and is the mass spurion defined 
by 



in terms of the light-quark mass matrix 



V 



and ^ — e''^/-'" . The covariant derivative is 

V^B = d^B+[V^,,B], 
and the vector and axial fields are 



(9) 
(10) 

(11) 

(12) 
(13) 
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The masses of the octet and decuplct baryons were first 
determined in Ref. [44] to NLO, 0(M|-_^). The octet 
baryon masses were later determined to NNLO, 0{M'^^^) 
in Refs. [45-47], and the decuplet baryons to NNLO in 
Refs. [48, 49]. 

The lattice calculation of baryon masses used for this 
work was a mixed-action calculation [27]. Therefore, 
the calculation and the corresponding low-energy ef- 
fective field theory for the baryons are inherently par- 
tially quenched. The partially quenched Lagrangian 
for heavy baryons was determined by Chen and Sav- 
age [50] through NLO. The extension of the partially 
quenched Lagrangian to the mixed-action Lagrangian, 
which includes corrections from the lattice spacing, also is 
known [51, 52], through NNLO [47, 49, 52]. At NLO, the 
corrections from the mixed-action are straightforward, 
amounting to corrections to the mixed valence-sea meson 
masses [53] appearing in the one-loop self-energy correc- 
tions. Fig. 1, as well as the partially quenched hairpin 
corrections [50]. A rigorous extrapolation to the con- 
tinuum limit requires multiple lattice spacings, but the 
lattice data used in this work was performed at a single 
lattice spacing (albeit with the size of the discretization 
corrections expected to be small, see Ref. [27] for details). 
However, one can make a qualitative estimate of the size 
of the discretization errors by comparing the extrapola- 
tions performed with the continuum extrapolation for- 
mulae and the mixed-action formulae. We now perform 
these chiral extrapolations. 

A. Variable projected minimization 

To construct the for the light-quark mass extrap- 
olation, we use the bootstrap resampled lattice data to 
form the covariance matrix, 

bs 

(14) 

where g is a supermass index which runs over both the 
mass relations M1-M7 as well as the different lattice en- 
sembles m007-m050.^ We then construct the x^, 

- E - /(M?' A)) C-^, (U,. - /(M,,, A)) , 

q,q' 

(15) 

where /(M^, A) is the HBxPT extrapolation formula for 
the mass relation Mg and depends upon the low-energy 



To the order we work in SC/(3) heavy baryon xPT, the relation 
Ms is exactly zero. In the partially quenched theory, this mass 
relation exactly vanishes at NLO, and at NNLO, there are resid- 
ual partially quenched effects which do not exactly cancel [49]. 



constants A = {Mq, Aq, bo, bp, bg, a, 7m, D, FjC^H,. . . }. 
At NLO, the baryon mass extrapolation formulae are lin- 
early dependent upon all the A,, except for D, F, C and 
7i, but they are linearly dependent upon and Ti.^. 
Therefore, to perform the numerical minimization, we 
first perform a linear least squares fit on all the linear 
LECs, \iin = {Mo, Ao,6D,feF,6o,CT,7M,C2,H^}, solv- 
ing for the low-energy constants as functions of D and 
F, and then perform the numerical minimization. This 
procedure is known as the variable projection (VarPro) 
method [54]. 

We first perform the minimization using relations 
M1-M7 and the NLO extrapolation formulae. In Ta- 
bles III and IV, we collect the results of both the contin- 
uum extrapolation and the mixed-action extrapolation. 
In these tables, the partially quenched/mixed- action low- 
energy constants a a/, /3 a/, ca/, 7a/ and cta/ are related 
to the LECs of SU{3) HBxPT of Eq. (8) by 

bD^\ (aM - 213m) , IM = 7if , (16) 

bp^l (5aA/ + 2/3a/) , = 4/"^^ > (17) 

1 2 

ho = (JM + -aM + -Pm ■ (18) 
6 3 

Then, for example, a a/ = um/Bq, where at leading order 
Alf^ = Bo^rris + mud)- In both cases, we perform the 
fit for four different ranges of the quark masses. The 
first fit includes only the lightest two mass points, while 
each successive fit includes an additional mass point. The 
results of the mixed-action and continuum fits are fairly 
consistent, as can be seen in Tables III and IV. 

Consistent with Refs. [27, 41], we find the values of D 
and F are significantly smaller than those determined ei- 
ther phenomenologically [42] or with the recent lattice de- 
termination [43] (these determinations provide D ~ 0.72 
and F ~ 0.45). With the use of the VarPro method, 
we have reduced the numerical minimization to a two- 
dimensional problem, thus allowing us to plot the result- 
ing x^ as a function of D and F . In Fig. 15, we provide 
a contour plot of the resulting x^ for the continuum fit 
using the lightest two ensembles. In Fig. 16, we plot the 
resulting x^ for the continuum extrapolation using the 
lightest five ensembles. Figures 17 and 18 display the 
same plots constructed with the mixed-action extrapola- 
tion formulae. In all these plots, the dark (blue) central 
area satisfies x^ ^ Xmin + 2.30, the ztcr confidence region 
for two parameters. Each successive n'^ contour repre- 
sents the ibncr confidence region. 

One observes that the mixed-action analysis returns 
values of D and F which are slightly larger than the con- 
tinuum SU{i) analysis, but the values are still smaller 
than expected. In Table V, we use the values of the 
LECs from the continuum analysis to compare with the 
experimental values of the individual baryon masses and 
the 1/Nc mass combinations in Table I (within a few 
sigma, the mixed-action extrapolations are consistent). 
Despite the apparent discrepancy in the values of D 
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TABLE III: Fit to Mass Relations Ri~Rr using NLO continuum HBxPT with a variation projection (VarPro) method. 



Fit Range 


D F 


Mo Ao &M pM 7m 5- a C'^ H'^ 
[GeV] [GeV] [GeV]-^ [GeV]-^ [GeV]"^ [GeV]-^ [GeV]-^ 


m007 - mOlO 
m007 - m020 
m007 - m030 
m007 - m040 


0.04(28) 0.14(6) 
0.17(8) 0.07(5) 
0.21(5) 0.08(3) 
0.20(5) 0.09(2) 


0.999(49) 0.21(16) -0.58(7) -0.42(16) 0.68(37) -0.10(6) -0.36(31) 0.01(12) 0.32(55) 
0.972(18) 0.48(5) -0.67(2) -0.81(5) 0.39(12) -0.11(2) 0.09(9) 0.21(3) 0.09(19) 
0.972(13) 0.44(5) -0.67(2) -0.76(3) 0.49(12) -0.12(1) -0.00(8) 0.15(2) 0.18(18) 
0.960(11) 0.43(04) -0.67(1) -0.73(3) 0.62(8) -0.13(1) -0.05(6) 0.14(2) 0.35(13) 



TABLE IV: Fit to Mass Relations R\~R-j using NLO mixed-action HBxPT with a variation projection (VarPro) method. 



Fit Range 


D F 


M) Ao ctM Pm 7a/ 5- a 
[GeV] [GeV] [GeV]"^ [GeV]-^ [GeV]-^ [GeV]'^ [GeV]"^ 


m007 - mOlO 
m007 - m020 
m007 - m030 
m007 - m040 


0.31(9) 0.26(5) 
0.36(3) 0.19(2) 
0.35(2) 0.16(2) 
0.34(2) 0.15(1) 


0.941(42) 0.242(73) -1.29(1) -0.46(1) 3.8(9) -0.35(5) -1.8(4) -0.043(9) 4.3(1.1) 
0.989(15) 0.365(30) -1.03(1) -0.70(1) 0.66(22) -0.24(2) -0.22(11) -0.009(4) 0.21(26) 
1.001(11) 0.368(28) -0.92(1) -0.71(1) 0.58(21) -0.20(1) -0.15(10) -0.007(4) 0.11(25) 
0.995(8) 0.376(26) -0.87(1) -0.70(1) 0.69(18) -0.20(1) -0.16(8) -0.001(2) 0.21(21) 



x' (D,F): m007 - mOlO NLO 5C/(3) 




FIG. 15: Contour plot of x^{D, F) constructed using a contin- 
uum 517(3) HBxPT x^ fit to the m007-m010 data sets. The 
dark (blue) inner region represents x^ ^ X™™ + 2.30, the 
±1(T confidence region for two fit parameters, D and F. Each 
successive n"^ contour represents the ±ncr confidence region. 



(D,F): m007 - m040 NLO 5(7(3) 




D 



FIG. 16: Contour plot of x^{D, F) constructed using a contin- 
uum 517(3) HBxPT x^ fit to the m007-m040 data sets. The 
dark (blue) inner region represents x^ ^ Xmin + 2.30, the 
±la confidence region for two fit parameters, D and F. Each 
successive n"^ contour represents the ±ncr confidence region. 



and F, one observes the extrapolated values of both the 
baryon masses and the mass combinations, measured as 
a percent deviation, are in reasonable agreement with 
experiment. 

We caution that one cannot draw strong conclusions 
about SU{3) HBxPT from this analysis, in particular 
the small values of D and F. First, the values of Mk,-k 
used in this work arc still larger than desirable for per- 
forming chiral extrapolations (in fact the strange quark 
mass is known to be ~ 25% too heavy [36]). Second, at 



NLO, the baryon masses are predicted to have large ^ 
corrections; however, this large nonanalytic contribution 
is not observed in the numerical results for the baryon 
masses themselves, see Fig. 2. Consequently, there must 
be strong cancellations between the different orders in the 
chiral expansion to produce the observed results. In or- 
der for the chiral extrapolation analysis to be consistent 
with both the baryon masses and the expected values of 
D and F ^ one may need to use lighter quark masses and 
to include 0{M^ ^) contributions to the masses. Fur- 
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^ (D,F): m007 - mOlO NLO MA 




D 



FIG. 17: Contour plot of x^{D, F) constructed using a mixed- 
action SU{6\3) HBxPT St to the m007-m010 data sets. 
The dark (blue) inner region represents ^ Xmi:n + 2.30, the 
±lcr confidence region for two fit parameters, D and F. Each 
successive n"^ contour represents the ±ncr confidence region. 



X' iD,F): m007 - m040 NLO MA 




D 



FIG. 18: Contour plot of x^{D, F) constructed using a mixed- 
action SU{&\2,) HBxPT fit to the m007-m040 data sets. 
The dark (blue) inner region represents ^ Xmin + 2.30, the 
±la confidence region for two fit parameters, D and F. Each 
successive n"^ contour represents the ±na confidence region. 



ther, it is likely that a combined analysis of both the 
masses and the axial couplings will be required. 

Performing the complete NNLO analysis of the octet 
and decuplct baryon masses introduces 19 new unknown 
LECs (for 30 total). Twelve of these LECs correspond 
to ©(to^) operators of the form Tt{BM+M+B) and 
give rise to K-ql corrections to the baryon masses. 
The other seven LECs correspond to interactions of the 
baryons with the axial current, of the form ^xiBA-AB)^ 
and give rise to baryon mass corrections of the form 
^^,/f,7,l'^(^l,K,r,)Mx- ^^^^ lattice data used in this 
work, Ref. [27], arc not sufficient to precisely determine 



TABLE V: Predicted baryon masses and mass relations from 
fit in Table III. The predicted values are listed according to 
the heaviest light-quark mass used. The m030 fit is consistent 
with that of m020 and m040. 





Exp. 


mOlO 


m020 


m040 




[GeV] 


[GeV] 


[GeV] 


[GeV] 


Miv 


0.939 


1.039(31) 


1.023(12) 


1.020(8) 


Ma 


1.116 


1.159(22) 


1.145(9) 


1.142(6) 


Me 


1.193 


1.219(22) 


1.231(9) 


1.221(7) 


Mh 


1.318 


1.306(16) 


1.309(7) 


1.303(5) 


Ma 


1.232 


1.376(30) 


1.454(12) 


1.427(9) 


Me. 


1.385 


1.461(28) 


1.531(10) 


1.516(8) 


Mh- 


1.533 


1.543(27) 


1.600(10) 


1.598(7) 


Mn 


1.672 


1.622(27) 


1.663(9) 


1.672(6) 


Ml 


175 


179(4) 


176(2) 


176(1) 


Ma 


-9.2 


-11(1) 


-13.6(4) 


-13.2(3) 


Ms 


-9.03 


-6.6(6) 


-7.5(2) 


-7.2(1) 


M4 


-0.21 


-0.15(1) 


-0.22(1) 


-0.20(1) 


Ms 


-0.21 


0.05(12) 


0.07(5) 


-0.00(4) 


Me 


-0.73 


-0.12(14) 


0.17(7) 


0.22(4) 


M7 


-0.092 


0.01(3) 


0.09(1) 


0.109(9) 


Ma 


-0.024 


0.00(1) 


0.009(5) 


0.001(4) 


Ms 


0.0097 


-0.003(8) 


-0.001(4) 


0.004(3) 


Mc 


0.0031 


-0.000(1) 


-0.0033(5) 


-0.0039(3) 



all these LECs. Further, the number of LECs for which 
the mass relations arc nonlincarly dependent increases to 
four. In performing the full analysis, we use the VarPro 
method to reduce the number of LECs minimized numer- 
ically to 9. To perform the partial NNLO analysis, we 
add only the O(m^) operators and are able to reduce the 
LECs which require numerical minimization to D and F ^ 
as in the NLO analysis. 

In Figs. 19 and 20, we display contours as a func- 
tion of D and F analyzed for the continuum S\J(?>) and 
the mixed-action extrapolation formulae, respectively. In 
both cases, the five lightest ensembles, m007-m040, are 
used. As can be seen, although the values of D and F 
arc still consistent with zero, they are now also consistent 
with the phenomenological and lattice values, albeit with 
large uncertainties. This large allowed variation of D and 
F translates into about an order of magnitude increase in 
the uncertainties of the predicted baryon mass relations 
as compared to those in Table V. The full NNLO analy- 
sis yields results in qualitative agreement with the partial 
NNLO minimization; however, there are more than two 
LECs upon which the x^ depends nonlincarly, and thus 
wc cannot display the corresponding contour plots. 

This analysis makes it seem plausible that a combined 
NNLO analysis of the baryon masses and the axial cou- 
plings could provide much more stringent constraints on 
the values of D and F ^ as well as the other LECs. Such an 
analysis would allow for the first rigorous exploration of 
octet and decuplct properties using S'C/(3) heavy baryon 
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{D,F): m007 - m040 NLO + NNLO c.t. 5(7(3) 




0.2 0.4 0.6 0.8 

D 



FIG. 19: Contour plot of x^{D,F) constructed using a con- 
tinuum SU(3) HBxPT fit including partial NNLO coun- 
terterms, as discussed in the text, to the m007-m040 data sets. 
The dark (blue) inner region represents ^ Xmin + 2.30, the 
±1(T confidence region for two fit parameters, D and F. Each 
successive -n}^ contour represents the ±na confidence region. 



x' (D,F): m007 - m040 NLO + NNLO c.t. MA 



^ 0.4 




FIG. 20: Contour plot of x^{D, F) constructed using a mixed- 
action SU{6\3) HBxPT x^ fit including partial NNLO coun- 
terterms, as discussed in the text, to the m007-m040 data sets. 



The dark (blue) inner region represents x 



Xr7 



, -f 2.30, the 



±lcr confidence region for two fit parameters, D and F. Each 
successive n**^ contour represents the incr confidence region. 



chiral perturbation theory. 



VI. CONCLUSIONS 

In this work, we have explored 1/Nc baryon mass re- 
lations using lattice QCD. Some of these relations were 
derived assuming approximate SU{3) flavor-symmetry, 
Mi-Mg, while other relations were derived assuming only 
SU{2) flavor-symmetry, Ma-Md- In all cases, we have 
found that the baryon mass relations are of the size ex- 
pected by both the 1/Nc power counting as well as the 



power counting in the SU(3) breaking parameter, e = 
(M|.-Af^)/A2 . For the mass relations Ma-Md, it would 
be interesting to consider strange quark values even heav- 
ier than those used in this work. Two of these relations 
become equivalent heavy-quark baryon mass relations as 
nis — > oo, and this exploration would probe the transition 
region between ms/Aqco < 1 — > 'tIj/Aqcd > 1. In the 
charm and bottom baryon spectra, these 1/Nc relations 
have led to stringent predictions for various mass com- 
binations. It also would be interesting to study Mi-Mg, 
which were derived assuming approximate SU (3) flavor- 
symmetry, as the strange quark mass is increased, to de- 
termine the value of nis at which the perturbative SU (3) 
expansion breaks down. 

Because of the definite spin and flavor transformation 
properties of the baryon mass relations, the discretization 
corrections to the mass relations also obey the 1 /Nc and 
e power counting, and arc suppressed. In fact, only Mi 
is subject to the leading 0{a^) mass corrections. The 
suppression of these discretization errors leads to more 
stringent constraints on the SU{3) heavy baryon chiral 
extrapolation of the mass relations to the physical point, 
and a better method for chiral extrapolation. It allows 
one to rigorously test the predictions and convergence of 
SU{3) HBxPT. 

In this work, wc also have performed both SU (3) and 
mixed-action chiral extrapolations of 1 /Nc mass relations 
using a variable projected minimization. Consistent 
with other works, including Ref. [27] (from which our nu- 
merical data are derived), we find that a NLO analysis, 
which includes the leading nonanalytic mass corrections, 
returns values of the axial couplings which are signifi- 
cantly smaller than expected from either phenomenology 
or lattice QCD. However, the predicted values of the mass 
relations, as well as the octet and decuplet masses them- 
selves, are in good agreement with experiment. Further, 
we have shown that the partial inclusion of NNLO mass 
corrections returns values of the axial couplings D and F 
which are consistent with expectations, albeit with large 
error bars. 

Unfortunately, the data set we have used in this work is 
insufficient to precisely constrain all the LECs in the full 
NNLO analysis. However, we have demonstrated that 
a simultaneous extrapolation of both the hypcron axial 
couplings as well as the 1 /Nc mass relations would allow 
for a rigorous exploration of the predictions and conver- 
gence of SU (3) heavy baryon xPT. It is possible, perhaps 
even likely, that the strange quark is too heavy for the 
SU{3) theory to be convergent. This outcome, however, 
fails to explain the successes of flavor SU{3) symmetry 
observed in nature in the baryon sector. Thus, a detailed 
study of this phenomena is warranted. 
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